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Abstract. 

The steady-state currents and densities of a one-dimensional totally 
asymmetric exclusion process (TASEP) with particles that occlude an integer 
number (d) of lattice sites are computed using various mean field approximations 
and Monte Carlo simulations. TASEP's featuring particles of arbitrary size are 
relevant for modeling systems such as mRNA translation, vesicle locomotion 
along microtubules, and protein sliding along DNA. We conjecture that the 
nonequilibrium steady-state properties separate into low density, high density, an 
maximal current phases similar to those of the standard (d = 1) TASEP. A simple 
mean field approximation for steady-state particle currents and densities is found 
to be inaccurate. However, we find local equilibrium particle distributions derived 
from a discrete Tonks gas partition function yield apparently exact currents within 
the maximal current phase. For the boundary-limited phases, the equilibrium 
Tonks gas distribution cannot be used to predict currents, phase boundaries, 
or the order of the phase transitions. However, we employ a refined mean 
field approach to find apparently exact expressions for the steady state currents, 
boundary densities, and phase diagrams of the d > 1 TASEP. Extensive Monte 
Carlo simulations are performed to support our analytic, mean field results. 
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1. Introduction 

There has been much recent interest in Asymmetric Exclusion Processes (ASEPs), 
with both periodic and open boundaries [1, 2, 3, 4]. In an exclusion process, particles 
move on a one dimensional lattice. With each movement a particle may move one 
lattice site to the right with rate p, or to the left with rate q. A particle may move 
to a lattice site only if that site is not already occupied by another particle. For a 
lattice with open boundaries, particles are injected into the chain with rate a only if 
the first (leftmost) site is empty. A particle occupying the last site is extracted with 
rate f3. Similarly, if reverse steps are allowed (partially asymmetric), injection into the 
rightmost lattice site occurs with rate <5 while extraction from the leftmost lattice site 
happens with rate 7. The steady-state currents, particle densities, and correlations 
of the ASEP have been studied extensively using both exact matrix product methods 
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and mean field approximations [1, 2, 5, 6]. The exact steady-state currents in the 
infinitely long lattice limit are summarized by a diagram depicting three regions; the 
maximal current, entry limited current, and exit limited current phases. Additional 
sub-phases corresponding to different density profiles, or domain walls, arise within 
the exit and entry limited phase regions [7, 8]. 

ASEPs have been used to gain insight into the noncquilibrium statistical 
mechanics arising in practical settings such as traffic flow [9, 10], ion channels [11, 12], 
mRNA translation [13, 14], and vesicle translocation along microtubules [15]. The 
standard ASEP assumes fixed transition rates for particles whose step sizes equal their 
diameters. However, in numerous applications, the driven particles and the underlying 
length scales of the microscopic transitions (steps) are not identical; i.e., the particle 
may be large and occlude d > 1 lattice sites. Specific biophysical examples include 
ribosomes moving along mRNA and large molecules or vesicles that are shuttled by 
motor proteins along microtubules. Ribosomes are ~ 20nm in diameter, but move by 
~ 2nm steps, codon by codon. Motor proteins that carry vesicles (> 50nm) typically 
move with ~ 5nm steps along microtubules. For these applications d « 10. 

Steady-state properties of the d = 1 ASEP (such as currents, densities, and their 
correlations) have been computed exactly using matrix product and recursion methods 
[5], When the size of the particles is greater than the size of a single lattice site, 
the matrix product method is inapplicable; however, diffusion constants and particle 
velocities have been calculated for two particles of arbitrary size [17]. Dynamical 
exponents have also been derived using Bcthe ansatz solutions [18]. In this paper, we 
use mean field analyses to derive apparently exact, simple formulae for steady state 
currents and densities of a totally asymmetric (q = 0) exclusion process (TASEP) 
with d > 1. Results for partially asymmetric processes, where particles can move in 
both directions, are straightforwardly derived from similar methods and will not be 
presented here. For the TASEP, we begin with the ansatz that the qualitative features 
of the current phase diagram remain intact when d > 1. Namely, that the steady-state 
currents exist in one of three regimes, high density (with currents determined by the 
exit rate), low density (with currents determined by the particle entrance rate), and 
the maximal current phase (with current determined by the interior hopping rate). 

In the following section, we consider a simple mean field theory that does not 
explicitly take into account the infinitely repulsive interactions of the large particles. 
The discrepancy between these qualitative results and those derived from Monte Carlo 
simulations are evident. We then derive quantitatively accurate analytical approaches, 
based on discrete Tonks gas partition functions and explicit state enumeration near the 
lattice boundaries. Results from both simple mean field, Tonks gas distribution, and 
the refined mean field approaches are compared with Monte-Carlo (MC) simulations 
in Section 3. We summarize and discuss extensions of our approach to partially 
asymmetric exclusions in Section 4. 

2. Mean Field Theories 

Consider identical particles of diameter d > 1 (in units of lattice sites) that are driven 
through a one-dimensional lattice having N 3> d sites. The particles interact through 
hard core repulsion. The ratios of entrance and exit rates, and steady state currents 
to the rate p are denoted a, (3, and J, respectively. With each motion, the particles 
move one lattice site to the right. 
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2.1. Simple Mean Field 

The simplest mean field approximation to the steady-states of this system assumes 
that the particles are on average uniformly distributed along the lattice without anti- 
correlations arising from their hard core repulsive interactions. The probability Pi 
that a particle is at site i is assumed to be independent of i. In this discussion, a 
particle's position in the lattice is determined by the location of its left edge. The 
left edge of a particle of size d can move forward only if the d sites ahead of it are 
unoccupied. Setting p equal to one (rescaling all rates in units of p), the steady state 
current is then: 

J = Pi(l - P i+1 )(l - P i+2 ) ■ • • (1 - P l+d ). (1) 

Away from boundaries, the locally uniform probabilities Pi ps p give 

J = p(l-p) d . (2) 

The maximum current that this uniform state can sustain, and its associated density 
is found by setting dJ/dp = 0, which yields 

1 d d 

Pmax = Jj-^Yy J max = (d + 1)d +l " ( 3 ) 

Now consider particle entry at the left end of the lattice, or particle exit from 
the right end. Figure 1 shows two extreme limits for the rules of entry and exit when 
d > 1. For example, a particle may only enter incrementally, one lattice site at a time 
(Fig. la), or it may enter completely in one step when the first d sites are empty 
(Fig. lb). One can also easily define an intermediate injection scheme appropriate for 
the specific application at hand. For example, in mRNA translation, the entry site 
for ribosomes must be completely cleared before another ribosome can bind. Similar 
microscopic considerations also permit different exit schemes (Figs. lc,d). However, 
in the N — ► oo limit, it is easy to see that the steady-state current and densities 
are independent of the different entrance/exit schemes. For example, partial entry 
is equivalent to complete entry into a lattice that is extended by d — 1 sites at the 
entrance boundary, while partial exit is equivalent to complete, immediate exit of a 
lattice extended by d — 1 sites at the right boundary. Since d/N — > 0, the differences 
arising from the variations in entrance and exit are finite-sized effects that vanish in 
the macroscopic limit. For the sake of simplicity and consistency, we will assume 
complete entry and exit (Figs. lb,d) in all of our subsequent work. 

For the entry-limited current, we set 0:^=1(1 — Pj) = Pi IljLi(l — Pi+j) an d 
approximate Pi+j ~ Pi- The occupation probability near the left entry site is thus 
Pi ^ a and the entry-limited steady-state current becomes 

J L = a(l - (4) 

which we obtain if ft is large enough (so that extraction is not the rate-limiting step) 
and a(l — a) d < i^wi , or equivalently, when a < l/(d + 1). Similarly, when the 

current is exit rate limited, /3P N -d+i Yl^Zl^-PN-d+i+j) = PN-dJl j= iO--PN-d+j)- 
For slowly varying occupation probabilities near the end of the chain, P/v ~ 1 — ft, 
and 

J R = ft d {\-l3). (5) 
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Figure 1. Entry and exit mechanisms, (a) Incremental entry, (b) Complete 
entry into the first d sites, (c) Incremental exit, (d) Complete, one-step exit as 
the particle reaches the last site. 



The extraction rate limited current (5) holds when (3 < d/(d + 1) and a > l/(d + 1). 
This simple mean field theory ignores the additional exclusion between particles with 
sizes greater than a single lattice site. The resulting steady-state currents and densities 
arc described by the following regimes: 
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The currents, boundary densities, and interior densities (6) reproduce the known 
results for the d = 1 TASEP. However, we do not expect these mean field results 
to produce accurate density profiles [2, 4]. Moreover, the boundary densities are 
expected to be inaccurate since they are independent of d. The predicted currents 
(6) are continuous across the a — (3 parameter space and the transitions are first 
order between phases (I) and (II) and second order between all other phases with a 
d-dependent discontinuity 
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(7) 



The currents, densities, and phase boundaries will now be compared with more 
accurate mean field theories and numerical simulations. 



2.2. Tonks Equilibrium Distribution 

We shall see that the simple mean field approximation developed in the previous 
section is inaccurate for d > 1. A more refined approximation for a TASEP with 
locally uniform particle distributions is to consider the Tonks gas in equilibrium. For 
example, in an infinite one-dimensional lattice at maximal current in the steady-state, 
the density distribution is uniform over the lattice except near the ends. 
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Figure 2. (a) Constraints and a configuration associated with the calculation of 
Z(n,N). (b) Typical configuration for the calculation of Z(n,N — 1). 



A fixed section of length N in the interior of the chain that contains n < N / d 
uniformly distributed particles can be treated as an equilibrium ensemble. This is 
because in an infinitely long chain, we can translate the segment TV at the mean 
particle velocity. This frame of reference allows us to treat the configurations sampled 
by the particle within N as that arising from the equilibrium grand canonical ensemble. 
In the large n limit, we can use the canonical partition function of a gas of hard rods 
on a discrete lattice (discrete Tonks gas) to compute the likelihood of configurations 
that support and contribute to the net particle current. 

Consider the case where the i th particle from the entrance of the lattice has its 
left edge located at site j (xi = j). The probability that lattice site j + d is empty (in 
other words, that the i th particle is free to move to the right) is calculated from the 
ratio P(xi + i > j + d\xi = j) = Z(n, N — 1)/Z(n 7 N), since Z(n, N — 1) corresponds 
to the number of states of n particles where particle i + 1 occupies sites to the right 
of j + d. One specific configuration in this sum is shown in Fig. 2a. 

The partition function for n particles of length d confined to a length of N > nd 
lattice sites is [191 



Z(n, N) 



N-(d- l)n 



(8) 



n 



which leads to 



P(x i+1 > j + d\x t = j) = 



Z{n,N - 1) 
Z{n, N) 

N - (d- l)n -n 1 - pd 



(9) 



(10) 



N-(d-l)n l- P {d-\y 
where p = n/N . The steady-state current in a uniform phase is thus 

J = P(xi = j)P{x l+1 > j + d\x t = j) = p 1 PC? 

1 -p(d-l) 

The maximal current J m ax and its associated density p ma x is found by setting d p J = 0, 
which yields 

Pmax = i— — 7= " and Jmax = " — 7= ~- (H) 
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Figure 3. Mean field scheme for particles of size d = 4 near boundaries, (a) d+ 1 
distinct states at the entrance region, (b) d + 1 distinct states at the exit region. 
All the states that are within d microscopic steps of the first state pictured are 
included. 



Equation (11) reduces to the standard TASEP result [1, 2] when d = 1 and 
Jmax, Pmax ~ l/<i as g? — > oo. The steady-state current will be the maximal current 
Jmax when the ASEP is forced under periodic boundary conditions (a ring) or as long 
as the injection and extraction rates of an open boundary ASEP are large enough for 
the internal moves to be the overall rate limiting steps. 

The Tonks gas approach would not be expected to yield accurate results in the 
boundary limited current regimes (I) and (II) since these currents are determined 
by the behavior at the boundaries where densities are not uniform. Even near the 
boundary limited and maximum current phase boundaries, when the density profile 
varies slowly (as power laws) [2, 4, 5], one cannot assume uniform density and apply 
a local (near the entrance and exits) Tonks gas. Nevertheless, within the maximal 
current phase, currents, and densities derived from the Tonks gas approach compare 
extremely well with those found from Monte Carlo simulations. 

2.3. Refined Mean Field Theory for Boundary-limited Currents 

Here, we develop a mean field theory that accurately models the behavior of the 
particles at the ends of the lattice, and matches with the Tonks gas theory as the 
injection a and extraction f3 rates arc varied across phase boundaries. The resulting 
steady-state currents are also found to agree, within error, with results from Monte 
Carlo simulations. 

First consider the entrance site at the left of the chain as shown in Fig. 3a. Upon 
balancing the probability currents into the first state shown in Fig. 3a, 

J L = aP{ Xl >d) = P{xx = d,x 2 > 2d). (12) 

Although the simple mean field assumption used to derive Jl and J a (Eqs. (4) and 
(5)) completely ignore exclusion beyond one lattice site, relative probabilities can 
be used to derive a refined mean field theory for the boundary limited currents. The 
exclusive nature of the particles are partially incorporated in this instance by explicitly 



7 



enumerating the d + 1 allowed states. For example, the probability that the first d 
sites are unoccupied is 



(1-0)° 

de{\ - ey- 1 + (i- 



p (xi >d)= Aal , md , (13) 



where is the relative likelihood that a site is occupied assuming a nearly uniform 
particle density. The normalisation factor d9(l — 9) d ~ 1 + (1 — 9) d is the sum of the 
weights of all possible states of the first d lattice sites (Fig. 3a). The probability that 
the first particle is at site d, and is not blocked from moving one lattice site to the 
right is: 

n^=d, X2 >2d ) = w -^f^ . (u) 

We have assumed a slowly varying density so that 9 represents a mean over the first 
d sites. The mean occupation of the first d sites is then 

d 0(1-0)*-! _ 9 

~! de(i - ey- 1 + {i - e) d i + (d-i)0 - 1 } 

Upon using Eq. (12), 9 = a, and 

a a(l — a) 

PL = — 7T— , Jl = : , ,, — jt- ■ (16) 
\ + {d—l)a 1 + (a — l)a 

For the exit rate-limited (small f3) steady-state current, we consider all possible 
states between sites N — d + 1 to N. Here, since the last particle lingers at site N, 
the occupation fractions will not be slowly varying over the last d sites (cf. Fig. 
7c). For small exit rates, we expect that close-packing of the last few particles 
near the end of the chain will split the densities into two branches. One branch 
corresponds to a lingering close-packed state with enhanced occupation fractions 9i 
at i = N, N — d, N — 2c?, .... These densities will be slowly varying allowing us to 
make the approximation 0jy ~ 9^-d ~ &N-2d ■ ■ ■ = @n- Sites belonging to the second 
density branch lie between the high occupancy sites. These second-branch sites tend 
to be occupied when the particles are not close-packed; for example, shortly after a 
particle is extracted from the lattice. These densities, 6m— i, 9n— 2, ■ • • , @N— d+i = 9 
are smaller than 9m by a factor of approximately f3. 

Transitions between states corresponding to occupancy in these two sets of sites 
are considered. Conservation of the probability of the first state (only site N occupied) 
involves the flux (3 out of that state, and the replenishment from the last state occurring 
at rate p. The steady-state current can be found from 

J R = 0P(x n = N) = P{ Xn = N — 1). (17) 

Since in each "unit cell" there is only one high occupancy type site, (site N), there is 
no additional normalization factor and the probability that the last particle is situated 
at site N, ready to exit with rate (3 is simply P(x n = N) = 9m- The probability that 
the last particle is at N — 1 is 



P{x -N-D- "(i-fl) d - 2 (i-M (18) 

[ " ' ~ (1 - 0)<<-i + (d - 1)9(1 - 9^-2 > 



8 



where (1 — 9) d 1 + (d — 1)0(1 — 0) d 2 is the total weight for all states of the sites 
N — 1, N — 2, . . . , N — d+ 1 (the last four rows in Fig. 3b). Upon solving (17), we find 

(—H = i (19) 

i-e N i + (d-2)e y ' 

To find another relationship between 9 and 9n, equate the flux into and out of the 
state where site N — 3 is occupied: 

Ml - M(l - Of 1 ' 1 = ^(l - #) d ~ 2 (l - M- (20) 
Upon solving equations (19) and (20), we find 

-,+ - *— ^ J W—Pl (211 

6n = pR i + {d-iw R i + (d-i)p' (21) 

and mean occupations at the low density sites pn-i = Pn-2, ■ ■ ■ ,PN-d+i = Pr = 
PPn = J R- 

The currents and densities (16) and (21) reduce to the exact results of the standard 
TASEP for d = 1. Unlike in the d = 1 case, particle-hole symmetry is no longer 
preserved for particles that occlude more than a single lattice site. A different state 
enumeration near the two ends of the lattice must be considered, but the resulting 
expressions for the currents Jl and Jr remain symmetric with interchange of a /3. 

2-4- Matching mean field phase boundaries 

The different approximations for the steady-state currents and densities cover all 
allowable values of the parameters a, (3. Equating the various currents (11), (16), and 
(21), we find the following phase boundaries between the different current regimes 

/y\ ^ * - 1 a T Q (l - oQ " 

(I) a<a= V3T? /J JL= l + (d-l)a pL= l + (d-l)a 

(II) ?<F=a*,a>0 Jn = P R =^0^ ^=1+1^ ^ 
(III) {a,p)>(a\p) J max = — 1— Pii 



(Vd+1) 2 ^ Vd(Vd+l)' 

Note that for d = 1 the boundary (and interior) densities predicted by all methods 
arc equivalent. However for d > 1 we will see that the current and boundary densities 
in equations (22) are significantly more accurate than those predicted by the simplest 
mean field theory (Eq. 6). In fact, the currents and boundary densities of (22) match 
our MC simulation results exactly. 

The qualitative difference between the steady-state boundary densities predicted 
by (6) and (16), (21) arise directly from the local, exclusion of the individual particles. 
The boundary densities are accurately predicted with mean field theory; however, we 
expect that the predicted density profiles are as inaccurate as they are for the standard 
d = 1 TASEP [4]. 

The phases (22) were constructed by compositing results from approximations 
unique to each of the three phases. The currents and densities are continuous at the 
phase boundaries. In fact, the predicted orders of the phase boundaries, as in the 
simple mean field theory, is consistent with the exact solution of the d = 1 TASEP 
[2, 4, 3]. Namely, that (dJ/da) and {dJ/df}) are continuous across the (I) - (III) and 
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(II) - (III) phase boundaries. The phase transition between (I) and (II) is first order 
while the transitions between (I) - (III) and (II) - (III) have the following d-dependent 
discontinuity in the second derivative: 

d 2 J 

da 2 „_„ 



as e 



G> 2 J 

da 2 



d 2 J 

d(3 2 



d 2 j 

dp 2 



2 



(23) 



The nature of the transitions are qualitatively different from those predicted using the 
simple mean field approach (7) and compare favorably with simulations. 



3. Monte Carlo Simulations 



Extensive Monte Carlo simulations were performed to validate the various analytical 
models presented in the previous two sections. As we expected the densities in the 
lattice to vary significantly as the injection and extraction rates were varied, we chose 
to base our Monte Carlo code on a continuous-time algorithm developed by Bortz, 
Kalos, and Lebowitz (BKL) [20]. The BKL algorithm has the advantage of maintaining 
a constant computational efficiency over a wide range of particle densities and has also 
recently enjoyed wide use in stochastic simulations of chemical reactions [21]. 

The finite-sized effects in our simulations was estimated by running with d = 1 
and lattices of varying lengths. For lattices at least one thousand sites long, the MC 
results were found to systematically deviate from the known TASEP results by less 
than half a percent. In the simulations with larger particles, the length of the lattice 
was scaled with d to ensure that at least one thousand particles could occupy the 
lattice simultaneously (N = 1000<i). The simulations were run for 4 x 10 9 steps, which 
was sufficient to reproduce the known TASEP results in lattices much longer than 
our d x 1000 site benchmark. In all simulations, a linear-congruential pseudorandom 
number generator with a period of 2 x 10 18 [22] was used. 



3.1. Currents 

Figure 4 plots the steady-state current as a function of particle size d in all three 
current regimes. While the simple mean field approximation always underestimates 
the currents, the Tonks distribution and the refined mean field (Eqs. (16) and (21)) 
give currents that match the Monte Carlo results within simulation error (« 1%) 

The currents for various fixed d as functions of a and (3 are shown in Figs. 5. 
From these curves one can estimate the critical values a* and (3* at which the phase 
boundaries lie. The mean field predictions for the currents displayed in Figs. 5 are 
a composite of the results from two different equations. Equation 11 was used for 
entrance and exit rates greater than the values of a* , (3* in Eq. (22). For smaller 
values of a,f3, Eqs. (16) and (21) were used. The simulated and refined mean field 
curves (22) are symmetric with respect to the interchange of a «-> (3. However, the 
simple mean field results (6) always underestimates the currents and are shown for 
comparison by the thin dotted, dashed, and long dashed curves corresponding to 
d = 2,4,8, respectively. Note the clear asymmetry between entrance and exit rate 
limited currents in the simple mean field approximation. 

Figures 4 and 5 demonstrate the extremely good agreement between the MC 
results and our optimal mean-field expressions (22) for the current. In all cases, the 
difference between the optimal mean field predictions and the MC results was within 
error. The extremely high quality of the fit between the mean field currents and the 




Figure 4. The currents within each of the phases (I), (II), and (III) as a function 
of integer particle size d. Results from simple mean field, refined mean field, and 
Monte Carlo simulations are shown, (a) The current Jj, as a function of d for 
f} = 10 and a = 0.1. (b) The maximal current J m ax when a = (3 = 10. (c). The 
exit limited current Jr for a = 10 and /3 = 0.1. The simple mean field approach 
always underestimates the currents since it smears out the density and disrupts 
the collective hop of each large particle. 




Figure 5. Currents as functions of a and f3 for fixed d = 2, 4, 8. (a) Fixed /3 = 10. 
(b) Fixed a = 10. The points correspond to MC simulations, the thick solid 
curves that fit the simulations represent the optimal mean field predictions (22), 
and the dotted, dashed, long dashed curves are the simple mean field predictions 
for d = 2, 4, 8, respectively. 
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Figure 6. Representative boundary densities in each of the three current regimes. 
< denotes the averaged (over the first d sites) linear density p\d at the left end 
when the current is in the entry rate limited (a = 0.1, fi = 10) low density phase 

(I) . • represents simulation results for the averaged density in the lattice interior 
under maximal current conditions (a = f3 = 1). ► represents simulation results 
for the density pjv at the very last site N. The thin continuous curves show the 
predicted densities from the Tonks gas and refined mean field approximations 

(II) , (16), and (21) and lie within simulation noise. 

MC results strongly supports our conjecture that the currents derived from the Tonks 
gas and refined mean field approaches are exact. 

3.2. Particle Densities 

We now compare the boundary densities predicted in (6) and (22) with those derived 
from MC simulations. The MC simulations also furnish density profiles which are 
difficult and tedious to compute using mean field approximations. In Figure 6, the 
left boundary and interior densities, multiplied by d (phd and PN/2d), and the right 
boundary density p~ft are shown as functions of particle size d. The net error in each 
data point was approximately 3%. The densities pLd, p^ i^d, and p\ from our refined 
mean field analyses (22) match the Monte Carlo results exceptionally well. 

We have also verified that all densities. pl,Pn/2> and P~r can be readily and 
accurately computed in all current regimes. For example, p^ can be found from 
aP(xi > d) = J, where J = J max for the maximal current regime (III), and J = Jr 
for the exit rate limited regime (II) . Similarly, the interior densities in the entrance or 
exit limited regimes are found by equating (10) to Jl or Jr in (22) and solving for 
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the appropriate root p. The boundary and interior densities for all regimes are 



(I) 
(II) 

(III) 



PL PN/2 Pr 

a i+(d-i)j L -y/(i+(d-i).j L )2-4dJ L a(l - a) /(3 
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d l + (d-l)/3 2d l + (d-l)/3 

(Vrf + i) 2 -i i i 
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Recall that in the standard d = 1 TASEP, mean field theory provides the correct 
boundary densities but incorrect density profiles [4]. Figure 7 depicts the full filling 
fractions pd as a function of site i for a TASEP with particle size d — 3. Notice a that 
a weak anti-correlation arises at the right boundary even in the entrance rate limited 
phase (inset of Fig. 7a). Similarly, there is a very small oscillatory behavior in the 
density near the entrance in the exit rate limited regime (first inset of Fig. 7c). A 
more striking anti-correlation is evident near the right boundary of the exit limited 
chain (second inset in Fig. 7c) where the density profile splits into approximately two 
branches with the branch corresponding to Pn, PN-d, PN-2d, ■ ■ ■ having the highest 
occupation. This is a consequence of the d-site exclusion strongly manifesting itself 
near the crowded exit and was used to obtain the refined mean field theory results (21). 
When currents are exit rate limited, the last particle in the chain spends a significant 
amount of time at the last site i = lOOOd. The next site available for occupation is 
at i = 999d thus forcing d — 1 sites to be empty as long as the last particle occupies 
the last site. Anti-correlations are negligible at the left boundary since all particles 
can only move to the right. In partially asymmetric exclusion models we expect the 
inclusion of backward motion would strengthen the anti-correlation effect at the left 
boundary of the lattice, and weaken it at the right. 



3.3. Phase Boundaries 

In order to map the phase diagram for the TASEP with finite-sized particles, the 
boundaries between the low density, high density, and maximal current phases were 
determined. Based on the good agreement between the MC currents and the currents 
predicted by equations (16) and (21), we chose to use these equations to predict 
the location of the critical values of a and (3. In the case of the boundary between 
the maximal current phase and the entrance/exit rate limited phases, this was 
accomplished by fitting the boundary limited current equations to the MC data. In all 
cases the fitting parameter was taken to be d, the size of the particles. The values of 
d produced by these fits were then used to find a*, (3* in eqs. (22). A variation of this 
fitting procedure was employed to determine the location of the boundary between 
entry and exit limited phases. Specifically, Jl from equation (16) with the fitted value 
of d was equated to a maximal current computed from the MC results. Equation (16) 
was then solved for a to produce an estimate for the critical value. 

To determine the order of the phase transitions, numerical derivatives of the 
current were computed using a fifth order centred difference formula. Although 
the results for the jump in the second derivatives d 2 J/da 2 and d 2 J/d(3 2 have noise 
magnified by the differentiation procedure, the results agree quite well (cf. Fig. 9) 
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Figure 7. Normalized density (pid) profiles in three representative current 
regimes, (a) Entrance rate limited regime (a = 0.1, /3 = 1). (b) Maximal current 
regime (a = (3 = 1). (c) Exit rate limited regime [a = 1, /3 = 0.1). The boundary 
regions are shown in the insets. 



with the predictions from the Tonks gas/refined mean field approach (23). An analysis 
similar to those that have been applied to domain wall dynamics [7, 23] is valid in 
the study of the TASEP with d > 1. For example, additional phase boundaries that 
delineate regimes of different domain wall dynamics can be found [7, 23]. One does not 
expect that larger (d > 1) particles to qualitatively change coarse-grained properties 
such as shock dynamics. Thus, the additional phase boundaries separating the regions 
of different density dynamics can be found by considering shock velocities of the form 
V = (J+ — J_)/(/9+ — p-) where J±,p± are the currents and densities to the right or 
left of a shock [7, 23]. 



4. Summary and Conclusions 



We have presented mean field theories which accurately model the currents and 
densities of a totally asymmetric exclusion process with particles of arbitrary size. 
Specifically, steady-state currents and densities were computed, using both the 
analytic, mean field approaches and Monte Carlo simulations. All calculations show 




Figure 8. (a) The phase diagram for the TASEP with particles of size d = 3 
derived from the Tonks distribution and the refined mean field (22). The circles 
represent the phase boundaries estimated from Monte Carlo simulations. For 
contrast, the d = 3 phase boundaries derived from the simple mean field model 
(6) are overlayed using thin curves and are asymmetric in a, (3. Boundaries 
across which first order phase transitions occur are shown with solid curves, 
while those across which second order transitions occur are drawn with dashed 
segments, (b) The critical values a* and /3* at which the (I) - (III) and (II) - (III) 
phase boundaries occur. For the Tonks gas and refined mean field approaches, 
a* = (3* = l/(\fd+ 1). Critical values derived from the simple mean field (6) are 
o* = l/(d + 1) and P* = d/(d + 1) > 1/2. 



that the steady-state currents qualitatively retain the three phase structure (low 
density high density, maximum current) familiar from the standard [d = 1) TASEP. 

Our best analytic approaches utilize an equilibrium Tonks gas distribution within 
the maximal current regime, and an exact enumeration of the most likely occupation 
states near the lattice ends in the entrance or exit rate limited regimes. Unlike the 
simplest mean field approach, these refined analytic approximations provide simple 
formulae for the steady-state currents and densities for all particle diameters d that 
agree extremely well (within simulation error) with our Monte Carlo simulation 
results. Based on this close agreement, we conjecture that the analytic expressions 
for the currents given by the Tonks/refined mean field approach (22) are exact. 
The phase diagrams produced by both of our analytic approaches and the Monte 
Carlo simulations are shown in Fig. 8. Both analytic approaches predict first 
order transitions between the high and low density phases and second order phase 
transitions between the maximal current and high/low density phases. However only 
the Tonks gas/refined mean field approach produces estimates for the discontinuities 
that match the simulation results. Although mean field methods are not expected to 
give correct density profiles near lattice boundaries, a consideration of the diffusion 
constant D = (1/2)(J + — J_)/(/? + — p_), of a coarse-grained density domain wall and 
the statistics of the domain wall position can be used to estimate the correlation or 
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Figure 9. The discontinuity in the second derivative of the steady-state current 
with respect to a or (3 at the maximal current phase boundary (I-III and II-III). 
The simple mean field approach (7) predicts an increasing discontinuity with d, 
while the Tonks/refined mean field approach (9) predicts a jump that decreases 
as 2/Vd- 



exponential density decay length £(<f) = \\n(J + (d) / J-(d))\ [7]. 

Since mean field results for d = 1 have been shown to be exact or partially 
asymmetric models [6, 24], we expect that our Tonks/refined mean field approach 
would also provide exact results for partially asymmetric exclusion processes with 
d > 1. One intriguing difference between d = 1 and d > 1 exclusion processes is the 
nonmonotonic behavior of steady-state densities near the exit. This behavior arises 
from the anti-correlations induced by particles of larger size and is a consequence 
of the total asymmetry in the dynamics. Although we have considered only total 
asymmetry, analogous expressions can be readily computed for partially asymmetric 
models by allowing backward particle steps (q ^ 0) and extraction/injection from the 
left/right. Mean field derivations for the partially asymmetric case would require the 
treatment of separate high and low occupancy sites near both ends of the lattice, as 
was done for the exit rate limited regime in the TASEP. 

The extension of exclusion processes to larger particles (d > 1) allows one to 
directly apply the them to many biophysical problems including ribosomc motion 
along mRNA and vesicle transport along microtubules. Our results, when generalized 
to include backward hops (q ^ 0), may be useful in determining the optimal scaling 
of step size and time increments in one-dimensional Brownian dynamics simulations. 
Furthermore, the Tonks gas distribution allows us to derive continuum dynamical 
equations near uniform states: 



• / s / \dp d 
P0M) = («-P)fl£fl£ 



P(l ~ pd) 
l-(d-l)p 



(25) 
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